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The impact of droplets on an inclined falling liquid film is studied experimentally using high¬ 
speed imaging. The falling film is created on a flat substrate with controllable thicknesses and 
flow rates. Droplets with different sizes and speeds are used to study the impact process under 
various Ohnesorge and Weber numbers, and film Reynolds numbers. A number of phenomena 
associated with droplet impact are identified and analysed, such as bouncing, partial coalescence, 
total coalescence, and splashing. The effects of droplet size, speed, as well the film flow rate are 
studied culminating in the generation of an impact regime map. The analysis of the lubrication 
force acted on the droplet via the gas layer shows that a higher flow rate in the liquid film produces 
a larger lubrication force, slows down the drainage process, and increases the probability of droplet 
bouncing. Our results demonstrate that the flowing film has a profound effect on the droplet impact 
process and associated phenomena, which are markedly more complex than those accompanying 
impact on initially quiescent films. 


I. INTRODUCTION 

Droplet impact on solid or liquid interfaces is a 
ubiquitous and fascinating phenomenon in nature [Dia, 
and has a wide range of applications; these include 
inkjet-printing, spray-painting, spray-cooling, internal 
combustion engines, fire suppression, deposition of solder 
bumps on printed circuit boards, surface-cleaning, and 
cell-printing. The outcome after a droplet impacts on 
a surface is a fascinating problem and has received 
considerable attention in the literature dating back to 
over a century ago [3]. Various impacting substrates 
have been employed, ranging from dry [4] and pre-wetted 
surfaces [5], surfaces with physical [6] and chemical 
heterogeneities mil], thin [9] and deep liquid layers [ini 
m, porous media [12], heated [nmu, sublimating m, 
and electrically-conducting [16], elastic membranes m, 
and granular materials m- Non-Newtonian [m Ho], 
surfactant-laden [21], and liquid-metallic droplets [2^1^ 
have also been studied. Different outcomes of impact 
have been observed, such as bouncing, coalescence, 
partial coalescence, splashing. The transitions among 
different regimes have been investigated 

Depending upon the choice of impacting surfaces and 
droplet physical properties, the impact dynamics can be 
extremely complex. Thus, analysing the accompanying 
phenomena during the impact process, and understand¬ 
ing the physics behind the phenomena are of critical 
importance. Generally, the impact dynamics of droplets 
is dominated by the interplay between surface tension, 
inertia, and viscous forces; this may be modified by 
the presence of non-Newtonian and/or physico-chemical 
and/or thermal effects. The relative importance of these 
forces can be quantified using dimensionless numbers, 
such as the Ohnesorge and Weber numbers. These di¬ 
mensionless numbers are used widely in the investigation 
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of droplet impact, e.g., in defining the regimes of impact 
[25] and in developing empirical correlations 0126]. 

Although droplet impact phenomena have been inves¬ 
tigated for more than a century, it is the recent devel¬ 
opment of high-speed imaging that allows researchers 
to unveil the fast evolving dynamic features of the 
phenomena [27]. Many phenomena occurring during 
impact remain poorly understood, and continue to 
attract interest. These phenomena include the collapse of 
the gas layer beneath the droplet during the impact [28]- 
|3Q] , the entrapment of air bubbles when a droplet impacts 
on a liquid surface im Eg, and the ejection of liquid 
sheets during the impact and its subsequent breakup 
into secondary droplets [33l [34]. Bird et al. [35] found 
that during the impact of droplets on angled or moving 
surfaces, the asymmetry leads to an azimuthal variation 
of the ejected rim, and the tangential component of 
impact can act to enhance or suppress a splash. Sikalo et 
al. [36] studied the impact of droplets on inclined dry and 
wetted surfaces, and found the bouncing on dry smooth 
or wetted surfaces, and no bouncing on rough surfaces. 
Gilet and Bush [37] studied the bouncing of droplets on 
an inclined substrate with a thin layer of high-viscosity 
fluid, and analysed the contact time and the conversion of 
energy. The liquid layer with high viscosity was used to 
ensure a smooth impact surface which facilitates droplet 
bouncing. Despite the large number of studies on droplet 
impact, the impact of droplets on flowing liquid films, 
to the best of our knowledge, has not been studied. 
This is surprising since it is a common phenomenon 
occurring in many two-phase flows, such as annular 
flows [38], spray-cooling [39], and spray-painting. The 
interaction between the droplet and the flowing liquid 
film could produce unique dynamic features distinct from 
the impact phenomena on quiescent liquids. This will be 
the subject of the present paper. 

A falling liquid film is used as the flowing liquid film for 
droplet impact in this study. Falling films are common 
in nature, important to many industrial applications, 
and are characterised by rich nonlinear dynamics gDl- 
EH] . The pioneering study of Kapitza was performed half 
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a century ago [44]. Since then, the flow in falling liquid 
films has been studied extensively through experimental 
measurements uni ES], low-dimensional modelling [ISl 
113 , and full numerical simulations [48]. The effects of 
many factors that influence the flow in falling films have 
been investigated, such as the effects of thermocapillarity 
[49] . electric fields [50] [51], and surfactants [52] [53] . 
Different processes that may be involved in falling films 
have also been studied, such as heat transfer [54], mass 
transfer chemical reactions [56], and phase change 
153 . Even though the flow of falling liquid films has 
been widely studied in the literature, it still remains an 
important research topic and attracts the attention of 
researchers of various disciplines. In this study, we use 
high-speed imaging to examine the impact of droplets 
of controlled diameter, angle of incidence and speed, on 
isothermal falling films in order to explore how the film 
dynamics influence the impact phenomena. In addition 
to the Ohnesorge and Weber numbers, this system is 
also parameterised by the film Reynolds number, which 
is based on the film mean thickness and velocity. The 
results of our work show that the film dynamics have 
a profound effect on the impact process, giving rise to 
phenomena that are not observed in the quiescent film 
case. 

The rest of this paper is organised as follows. The 
experimental system is described in Section [TT| including 
the falling film rig, the droplet generation unit, and the 
high-speed imaging system. The experimental results are 
presented in Section [nil ^ regime map of droplet impact 
is created, and the effects of the Weber, Ohnesorge, and 
Reynolds numbers are studied by varying the impact 
speeds, the droplet sizes, and the flow rates of the falling 
film. A simple model based on lubrication theory is 
used to explain the bouncing phenomenon. Finally, 
concluding remarks are provided in Section [TV| 


II. EXPERIMENTAL METHOD 
II. 1. Falling liquid film 

The experimental setup used to study droplet impact 
on flowing liquid films is shown in Figure The falling 
film rig comprised a rigid and impermeable plate 41 cm 
long and 28 cm wide made of a titanium foil with a 50 
/im thickness inclined at 45° to the horizontal. The liquid 
used was deionised water of density 10 ^ kg/m^, viscosity 
10“^ Pa-s, and surface tension 72.8 mN/m, and the flow 
rate was controlled by a valve and a variac transformer 
for the pump. An ultrasonic flow metre (Flownetix 100 
series, Flownetix, UK) was used to measure the flow 
rate of the falling film, which was in the range 0.9T0.I- 
I0.7T0.3 1/min. The lower limit of the flow rate of the 
liquid film is the minimum flow rate that can ensure the 
whole substrate is wetted. The mean film thicknesses 
h and the mean flow speed u were estimated using the 
Nusselt solution [43] to be in the range 0.29-0.65 mm, 
and 0.19-0.98 m/s, respectively. The Reynolds numbers 



FIG. 1. Schematic diagram of the experimental setup for 
droplet impact on falling liquid films which comprises a falling 
film rig, a droplet-generation unit, and a high-speed imaging 
system. 


of the falling films Re = pq/{wfi) were therefore in the 
range 54-636, where p and p are the density and the 
viscosity of the liquid, q and w are the flow rate and 
the width of the falling liquid film, respectively. Since 
laminar-turbulence transition in falling films happens at 
Re > 1000 [40], the flow in the falling liquid film studied 
here is laminar. 


II.2. Droplet generation unit 

Droplets of uniform size were generated from the tip of 
blunt syringe needles. Deionised (DI) water was pumped 
into the needle using a syringe pump (Braintree Scientific 
Inc., UK) at constant flow rates. The droplets broke up 
at the tip of the syringe needle in the dripping mode, fell 
under their own weights, and impacted obliquely onto 
the inclined liquid film. The dripping mode and the 
gravity-driven acceleration of the droplets ensure that 
the process is highly reproducible. To study the effect 
of droplet size on the impact process, the droplet size 
was varied via a change in the diameter of the needle. 
Three droplet diameters were studied, 3.3, 4.0, and 5.2 
mm with standard deviation less than O.I mm, while the 
corresponding Ohnesorge number,Oh = were 

0 . 0021 , 0.0019, and 0.0016, respectively. Here, p is the 
density of the liquid, 7 is the surface tension, p is the 
dynamic viscosity, and d is the droplet diameter before 
impact. To study the effect of the droplet speed on the 
impact process, the impact speeds of the droplets were 
varied by changing the vertical distance of the needle 
towards the falling liquid film from I to 50 cm. The 
corresponding speeds of the droplet vo were in the range 
0.30T0.02-2.99T0.09 m/s, and the corresponding Weber 
numbers. We = pdv^/^^ were in the range of 4-630 for 
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different droplet sizes. The impact point on the liquid 
film is 34 cm from the inlet of the falling film, which is 
larger than 500 times of the film thickness of the Nusselt 
solution m at the highest flow rate studied, and allows 
the development of the film flow on the substrate. The 
details of the flow in the falling film will be discussed in 
Subsection IIII.ll 


II.3. high-speed imaging system 

The droplet impact process was recorded using a high¬ 
speed camera (Olympus i-SPEED 3) at a frame rate 
of 5000 frames per second. A macro lens (Sigma EX 
105mm) was used to capture the details during the 
impact. Two halogen lamps (500W) were used for 
illumination. To minimize the heating effect of the 
halogen lamps, the lamps were kept on only when the 
camera was recording. Side-view images were captured 
at a resolution of 804x396 pixels (about 40 /im/pixel) 
with 10° deviation from the plane of the film, while top 
view images were capture at a resolution of 528x600 
pixels perpendicular to the film (about 40 /im/pixel). 
The high-speed images were analysed using customised 
Matlab routines. The speed and the size of the droplets 
were measured from the high-speed images. 


III. RESULTS AND DISCUSSION 

In this section, we present a discussion of our results 
starting with a brief description of the wavy falling film 
dynamics in the absence of droplet impact; this provides 
an indication of the expected state of the interface 
prior to droplet impact at various film Re. A detailed 
investigation of the wavy flow in the liquid film is beyond 
the scope of this study. We then describe the main 
phenomena associated with the impact process, before 
presenting a flow regime diagram as a function of system 
parameters. The results of a simple model based on 
lubrication theory are also discussed, and used to explain 
the trends associated with a subset of the phenomena 
observed. 


III.l. Wavy flow in liquid films 


Topview 





Flow 

direction 


5 mm 



FIG. 2. Typical images of the wavy flow in the falling liquid 
films at different Reynolds numbers: (a) Re = 54, (b) Re = 
245, and (c) Re = 636. The flow rates are 0.9, 4.1, and 10.7 
1/min, respectively. 


of consecutive waves, and the (temporal) frequency of 
the waves is the reciprocal of the horizontal distance of 
consecutive waves. The amplitude of the waves could not 
be measured in this experiment quantitatively due to the 
limitation of the experimental facility. The variation of 
the film thickness, however, is represented by the light 
intensity. Previous studies of falling liquid films [40H43] 
have shown that the amplitude of the waves is usually 
smaller or comparable to the mean film thickness, and 
the wave lengths are much larger than the mean film 
thickness. The surface curvature of the liquid film is very 
small as compared to the droplet curvature due to the 
small aspect ratio of the liquid film (thickness vs film 
length). At low Re, the solitary wave in Eigurej^ travels 
downstream at a constant speed. With increasing Re, 
the waves interact with adjacent waves as they propagate 
downstream, as shown in Eigure[^-c. The speed of the 
waves can be estimated from the gradient of the line of 
the waves in the time-space plot, which shows that the 
wave speed increases with Re. 

In the following subsections, the impact processes on 
falling films at different conditions of We, Re, and Oh 
numbers are analysed. Various combinations of the 
parameters result in different phenomena as observed 
in the experiments: bouncing, partial coalescence, total 
coalescence, and splashing, as will be detailed below. 
Because the waves on the liquid films are formed 
naturally and stochastically, the drop may impact on 
different parts of the waves. Therefore, the outcome of 
the impact may vary, and the results are presented in a 
probability manner. 


The wavy flow in the falling liquid films was analysed 
from the top-view and side-view images near the point 
of the droplet impact, as shown in Eigure At low 
Re values, the liquid films are characterised by interfaces 
that exhibit mild deformations with occasional solitary 
waves. These waves are formed with coherent precursor 
waves, as shown in Eigure [^. As Re increases, the waves 
become more frequent and less coherent. The evolution 
of the waves can be seen from the time-space plot of a 
line along the flow direction, which is obtained from the 
top-view images, as shown in Eigure In a time-space 
plot, the (spatial) wave length is the vertical distance 


III.2. Bouncing 

‘Bouncing’ refers to situations wherein a droplet 
approaching the falling film, does not coalesce with it, 
but instead bounces back, as shown in Eigure which 
was generated for Oh = 0.0021, We = 4, and Re = 636. 
The coalescence of a droplet with a liquid film involves 
three stages: (i) droplet approach towards the film; (ii) 
the thinning of the thin, intervening gas layer between 
the droplet and the film; and (hi) the rupture of the 
gas layer. The first stage is driven by the momentum 
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FIG. 3. Time-space plots showing the evolution of the wavy 
flow in the falling liquid films at different Reynolds numbers: 
(a) Re = 54, (b) Re = 245, and (c) Re = 636. The flow rates 
are 0.9, 4.1, and 10.7, 1/min respectively. The speeds of the 
waves can be estimated from the gradient of the dashed lines 
of the waves. 


of the droplet gained from gravity, the second stage is 
resisted by the viscous force in the thin gas layer, and 
the third stage is dominated by capillary force. As the 
gas layer becomes thinner, the resistance becomes larger, 
according to lubrication theory [58l[59]. Therefore, the 
gas layer could potentially stabilise the droplet against 
coalescence. Rupture of the gas layer occurs when its 
thickness decreases below a certain level so that van der 
Waals forces come into operation to destabilise the layer 

m- 

Prior to impact, the droplet is almost spherical (see 
Figure]^), as the capillary force minimizes the surface 
area and, in turn, the surface energy. Upon impact, the 
droplet progressively adopts a pyramidal (see Figure]^) 
and then a pancake shape (see Figure]^). The latter 
increases the contact area of the droplet with the film and 
further retards the thinning of the gas layer, according 
to lubrication theory [SI Eg. The surface energy of the 
droplet, can be estimated through image analysis 
by assuming a surface of revolution, = S'y, where 
S is the surface area of the droplet. The result shows 
that from the undistorted droplet shape (see Figure 
c) to the maximum deformation (see Figure 4g), the 
surface energy increases from approximately 2.5x 10“^ J 
to 3.4 X 10“^ J. Meanwhile, the kinetic energy before 
the impact, Ey = and the loss of potential 

energy. Eg = mgAz^ are approximately 8 x 10“^ J 
and 5 x 10“^ J, respectively. Here, m is the mass of 
the droplet, and Az is the vertical distance that the 
droplet travels from the instant of apparent contact to the 
instant of maximum deformation. Therefore, the relative 
conversion of energy from kinetic and potential energy to 
surface energy during the process can be calculated from 


AE^ 

E^+ Eg' 


( 1 ) 


which is estimated to be about 70%. Such high ratio 
of energy conversion has also been found for droplet 
bouncing on horizontal and inclined wet surfaces m 
The Weber number is a key parameter that affects the 
thinning of the gas layer prior to possible rupture: lower 
Weber numbers (lower impact speeds) lead to longer 
gas layer drainage times, and, therefore, an increased 
probability of droplet bouncing. The momentum of the 
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FIG. 4. ‘Bouncing’ of a droplet on a falling liquid him. The 
dimensionless parameters are Oh = 0.0021, We = 4, and 
Re = 636. The droplet speed is 0.3 m/s, the how rate of the 
falling him is 10.7 1/min, and the droplet diameter is 3.3 mm. 
See Ref. [62], Movie 1. 


falling him, on the other hand, for a hxed inclination 
angle, increases with the him Reynolds number (the 
liquid how rate). Therefore, droplets have a greater 
chance of bouncing on liquid hlms with higher Reynolds 
numbers. This is consistent with the phenomenon of 
droplet levitation at the hydraulic jump point of a 
liquid him produced from a vertical jet impinging on a 
horizontal surface m- The details will be explained in 
detail in Section IlII. 71 


III.3. Partial coalescence 

For a range of system parameters, the droplets were 
observed to undergo ‘partial coalescence’ with the falling 
him. This corresponds to situations wherein, upon 
impacting the him, a fraction of the droplet coalesces 
with the him, but immediately after impact, a small 
daughter droplet or a cascade of daughter droplets 
j63l l64) are created; an example of partial coalescence 
is shown in Figure It occurs as a consequence of small 
Weber numbers and small him Reynolds numbers. 

The rupture of the gas layer starts with the formation 
of a liquid bridge connecting the droplet and the liquid 
him. Due to the large curvature at the liquid bridge 
and the corresponding large capillary force, the liquid 
bridge immediately grows, and produces capillary waves 
on the surface of the droplet and on the liquid him 
(see Figure and Figure for the formation and the 
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FIG. 5. ‘Partial coalescence’ of a droplet on a falling liquid 
film. The dimensionless parameters are Oh = 0.0021, We = 4, 
and Re = 54. The droplet speed is 0.3 m/s, the flow rate of 
the falling film is 0.9 1/min, and the droplet diameter is 3.3 
mm. See Ref. Movie 2. 


propagation of the capillary waves). The waves on the 
droplet then propagate rapidly towards the apex of the 
droplets (see Figure [^). Immediately after the collapse 
of the gas layer, the liquid at the bottom of the droplet 
is significantly affected by the capillary waves and is 
pulled towards the liquid film. However, the liquid at 
the droplet apex is less affected, and there is no sufficient 
time to catch up with that near the bottom of the droplet. 
Therefore, the droplet becomes a long thin liquid thread 
(see Fignre[5l) , which eventually breaks up into two parts 
(see Figure!^) due to the inward momentum of the 
collapsing neck [63] . The lower part mixes with the liquid 
in the falling film, while the upper part gradually falls 
on the liquid film due to gravity as a daughter droplet. 
The daughter droplet tends to be more stable against 
coalescence because of its smaller speed than the primary 
droplet before partial coalescence, as shown in Figure [^. 

The surface energy after pinch-off of the daughter 
droplet (see Figure 5n) can be estimated through image 
analysis to be ^ 6x 10“^ J. It is much smaller than the 
initial surface energy before the impact (see Figure 5a), 
F^^^2.5xlO“^J. This indicates significant energy loss 
due to viscous dissipation and the formation of capillary 
waves during the process. 
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FIG. 6. Formation and propagation of capillary waves 
immediately after the collapse of the gas layer between the 
droplet and the liquid film. The dimensionless parameters 
are Oh = 0.0021, We = 4, and Re = 54. The droplet speed is 
0.3 m/s, the flow rate of the falling film is 0.9 1/min, and the 
droplet diameter is 3.3 mm. 



FIG. 7. ‘Total coalescence’ of a droplet on a falling film. The 
dimensionless parameters are Oh = 0.0021, We = 76, and 
Re = 636. The droplet speed is 1.3 m/s, the flow rate of the 
falling film is 10.7 1/min, and the droplet diameter is 3.3 mm. 
See Ref. [62], Movie 3. 


cence is shown in Figure [^ generated for Oh = 0.0021, 
We = 76, and Re = 636. If the Weber number is above 
a certain level, the upper part of the droplet, follows the 
lower part of the droplet and coalesces with the liquid 
film (see Figure [^), without breaking up into a daughter 
droplet, as in the case of partial coalescence. In addition, 
the high-momentum liquid in the droplet pushes away 
the liquid in the flowing film and produces a crater (see 
Figure]^). The waves produced by this event propagate 
away from the crater and dissipate gradually (see Figure 

Si)- 


III.4. Total coalescence 

In contrast to the situations considered above, ‘total 
coalescence’ refers to complete coalescence of the droplet 
with the falling film upon impact without forming 
smaller, secondary droplets; an example of total coales- 


III.5. Splashing 

If the momentum of the droplet is sufficiently large, 
characterised by high We, the droplet drives the liquid in 
the film to form a ‘crown’, which then disintegrates into 
many droplets, as shown in Figure [^ This ‘splashing’ 
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FIG. 8. ‘Splashing’ during droplet impact on a falling film. 
The dimensionless parameters are Oh = 0.0021, We = 316, 
and Re = 636. The droplet speed is 2.67 m/s, the flow rate 
of the falling film is 10.7 1/min, and the droplet diameter is 
3.3 mm. See Ref. [62], Movie 4. 


phenomenon occurs because the gas layer separating the 
impacting droplet and the falling film undergoes rapid 
rupture, and capillary forces are not sufficiently strong 
so as to balance the large force per unit area at high 
Weber numbers. Thus, following impact, the liquid 
undergoes severe deformation leading to the formation 
of a cylindrical sheet of liquid around the droplet (see 
Figure]^), which then assumes a crown shape (see Figure 
[^). The crown rim then undergoes an Rayleigh-Plateau 
instability [65] that results in the creation and subsequent 
ejection of small droplets (see Figure [^). The crown 
formation and the ejection of small droplets for the 
impact of droplets on quiescent liquids has been well 
studied in the literature [mm EH El]. Being different 
from impact on quiescent liquids, the wavy interface of 
the liquid film results in a more irregular crown structure 
and different sizes of the small droplets. The crown is also 
not symmetric due to the oblique impact angle [35]: the 
rim on the downstream side is much higher than that 
on the upstream side of the crown, which leads to the 
ejection of more droplets on the downstream side. The 
details will be discussed in Subsection IIII.9I 



FIG. 9. Regime map for droplet impact on falling liquid 
films. Here, ‘S’, ‘T’, and ‘M’ designate ‘splashing’, ‘total 
coalescence’, and multiple phenomena as described in Figure 
[To] The Ohnesorge number is fixed at Oh = 0.0021. The 
droplet speed is in the range of 0.30-2.99 m/s, the flow rate 
of the falling film is in the range of 0.9-10.7 1/min, and the 
droplet diameter is 3.3 mm. The lines are used only to guide 
readers’ eyes. 


III.6. Regime map 


A regime map for droplet impact on falling liquid films 
is created by varying the flow rate of the falling film and 
the impact speed of the droplets; this is presented in 
Figure as a Re vs We plot which shows the regime 
transitions as these parameters are varied for Oh = 
0.0021. For a given Re, large We values lead to splashing, 
while moderate We values lead to total coalescence. 
At the lowest We studied, multiple phenomena are 
observed, namely bouncing, partial coalescence, and 
total coalescence, and their probability of occurrence in 
repeated experiments are plotted in Figure [Toj 


As shown in Figure partial coalescence always 
occurs for the lowest Re studied. As Re increases, 
the probability of partial coalescence decreases, while 
that of total coalescence rises. This is because a thick 
liquid film facilitates the flow of liquid in the direction 
perpendicular to the film and allows the liquid near 
the top region of the droplet to catch up with that 
near the bottom region, thereby reducing the probability 
of formation of secondary droplets following impact. 
The probability of bouncing also increases with the 
film Reynolds number. This is consistent with the 
phenomenon of droplet levitation at the hydraulic jump 
point of a liquid film produced from a vertical jet 
impinging on a horizontal surface m, which corresponds 
to a large Re and negligible We in Figure 
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Re 

FIG. 10. Probability of occurrence of different phenomena as 
a function of Re with Oh = 0.0021 and We = 4. The droplet 
speed is 0.30 m/s, the flow rate of the falling film is in the 
range of 0.9-10.7 1/min, and the droplet diameter is 3.3 mm. 


film-gas interface in the present, approximate analysis: 


d'^Ug 

II 

(2) 

d^ui 

_ dpi 

(3) 


dx ’ 

d'^Ud 

_ dpd 

(4) 


dx ’ 


where the subscript g indicates the gas layer, and I the 
liquid film, and d the droplet. Then the velocity profiles 
along the 2 : direction are parabolic: 

Ug = AgZ^ + “ 1 “ ( 5 ) 


III.7. Lubrication gas layer and force analysis 


To explain the phenomena presented above, the gas 
layer between the droplet and the liquid film and the 
forces involved are analysed. Prior to the rupture of 
the gas layer between the droplet and the liquid film, 
the gas layer drains gradually. The drainage of the gas 
layer is controlled by the competition between the inertia 
of the approaching droplet and the lubrication force in 
the gas layer. The former serves as the driving force 
and is indicated by the normal component of the droplet 
speed to the falling film, while the latter serves as the 
resistance and is affected by the flow speed in the film. 
Thus, droplet bouncing is expected only at low normal 
speeds of impact. Even though the thickness of the gas 
layer is too thin to be determined from the observation in 
Figure measurements of the film thickness for droplet 
impact on solid surfaces in the literature [28l |66] have 
shown that the thickness of the gas layer can be down 
to several micrometres. To estimate the lubrication force 
acted on the droplet by the gas layer during the impact 
process, two complementary methods are presented. The 
first method is through an analysis of the flow in the gas 
layer, the liquid film, and the droplet, while the second 
method relies on the analysis of the droplet deformation 
shown on the high-speed images. 


III. 7.1. Estimation from the flow in the gas layer, the 
liquid film, and the droplet 

Due to the small thickness of the gas layer and the 
liquid film relative to the drainage length of the gas 
layer in the tangential direction, the flow in the gas 
layer and in the liquid film can be approximated using 
lubrication theory. The approximation is also applied in 
a penetration depth hd in the droplet, and the velocity 
gradient vanishes and the velocity approaches the droplet 
tangential speed at the penetration depth EH EH] , as 
shown in Figure We also neglect deformations of the 


— Aiz‘^ Biz Cl, (6) 

— AdZ BdZ + Cd-, (7) 

where the coefficients Ag, Bg, Cg, Ai, Bi, Ci, Ad, Bd, 
Cd can be obtained by satisfying the following boundary 
conditions and closure equations. A no-slip boundary 
condition is imposed at the substrate, and continuity of 
velocity is imposed at the gas/liquid-film interface and 
at the gas/droplet interface. 


«iL=o = 0’ (8) 

^9\z=h, = (9) 

^a\z=ho ~ (^ 9 ) 

where ho = hi hg. The matching of shear stress at the 
gas/liquid-film interface and at the gas/droplet interface 
yields 


hg 


dz 


= hi 


z=hi 


dui 

dz 


z=hi 


( 11 ) 


Ms 


dUg 

dz 


= hi 


z=ho 


dud 

dz 


z=ho 


( 12 ) 


The shear stress at the penetration depth of the droplet 
vanishes, and the speed approaches the tangential speed 
of the droplet with respect to the substrate, I/q- 


Ud\^=h2 = ^ 0 ’ (^ 9 ) 


dud 


= 0 , 

Z = h2 


(14) 


where /12 = hi hg hd- The flow rate of the liquid film 
is given by 


w 



uidz = q. 


(15) 












The pressure gradients in the liquid film and in the 
gas layer are assumed to be the same because film 
deformations are neglected, and normal viscous stresses 
correspond to high-order terms within lubrication theory; 


thus, Eqs. <§ and (§ yield 

d^Ug d‘^ui 


( 16 ) 


Therefore, the coefficients in Eqs. (i-® can be ob¬ 
tained: 


_ 3^iiD2 

^9 ~ 7 2 ’ 

0 = 0 , 

where 

Di = Sho^iJg - Sho^iJi - hi'^iJig + hi'^iJi - ‘ihoh2Hg - hohijig + h2hiHg + 2hohiiJi 

£>2 = {hoiig - h 2 Hg - 2hoiJ,i - 2hm,g + 2hiiii) q/w + Uohi'^i^g 

Ds = (3/io^Mg - 3/io^W - 3/iz^M5 + 3hi‘^i2i - 3/io/i2Mg) q/w + Uohi^Hg 

£>4 = ( 6 /io^W - 6 /io^Mg + 6/io/i2Mg + ihohi^ig - 3h2hifig - Ghohi/^i) q/w + hi^iigllo 
£>5 = {Sho^iJi - ^hi^jJLg + + QhohiiJg - QhohiiJi) q/w + {hi^/jg - Shohi'^jig) Uq 


_ 2^iiDz _ 

-^g — 2 r-k ’ ^9 — 


Ad = 


hi HgDl 

D5 


hi {ho — / 12 ) -Di 


Bd = 


{hg hi) B 4 

hlflgDl 

—2Adh2, 


Ai=- 


3 D 2 

hi^Dd 


Bi = 


2£>3 

hi'^Dd 


Cd — Adh\ + U 01 


(17) 


Although the exact value of pressure could not be known 
due to the complex structure at the rim of the gas layer, 
the pressure difference between the gas layer and the 
ambient can be estimated from the pressure gradient, 

ddpg d d'^Ug 
“ 2 ^ = “ 2 ^®^ = 

where we have made use of Eqs. ([^, (|^, ( p!^ , and ( [l7| ). 
Therefore, the lubrication force can be estimated as 

-^lubrication — {Pg Pa) dA ^ — d flgAg. (f9) 

A 

The liquid film thickness hi can be estimated from the 
Nusselt solution for the falling liquid film [43], 

, 20 ) 

\wpg cosU) 

The penetration depth in the droplet is assumed to 
be hd = as in Ref. [ 68 ]. Therefore, during 

the drainage of the gas layer, the magnitude of the 
lubrication force Tiubrication can be estimated at different 
thicknesses of the gas layer, as shown in Eignre[T^ The 
lubrication force (i.e., the resistant force) being large 
than the normal component of the gravitational force 
(i.e., the driving force) is a minimum requirement of 
droplet bouncing. Therefore, the gravitational force of 
the droplet in the direction normal to the falling film, 
-^gravity = COS (d), is uscd as a reference to scale 

the estimated lubrication force, where 0 is the inclination 
angle of the falling film to the horizontal. Eignre [T^ shows 
that as the gas layer becomes thinner, the lubrication 



FIG. 11. Schematic diagram of the flow in the gas layer 
and the liquid film beneath the droplet during the impact 
on flowing liquid film. 


force becomes larger. It also shows that the lubrication 
force increases with increasing the film flow rate. Since 
a large lubrication force indicates a large resistance of 
the droplet, a slow drainage of the gas layer, and a high 
chance of droplet bouncing, the result in Eigure is 
consistent with the experimental observation shown in 
Eigure [To] 


III. 7.2. Estimation from the droplet deformation in 
high-speed images 

The force applied on the droplet by the falling liquid 
film via the gas layer can also be estimated from the 
droplet deformation in the high-speed images. When 
the droplet reaches the maximum deformation, it has a 
pancake shape, as shown in Eigure [^. The pressure 
in the droplet can be estimated through the surface 
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FIG. 12. Estimation of the lubrication force normalised 
by the gravitational force normal to the him, Egravity = 
pg cos (0), as a function of the how rate of the liquid 
him at different thickness of the gas layer. 



P, 

Liquid film 


FIG. 13. Schematic diagram of the droplet deformation and 
the pressure. 

and viscous forces. The inertial force tends to deform 
the interface, while the capillary and viscous forces tend 
to resist this deformation. Their relative importance can 
be quantihed by their ratios represented by the Weber 
number and the Reynolds number for him deformation, 
respectively given by: 

(24) 

7 


curvature of the droplet. 


Pd 


1 1 
R^r 


+ Pa 


( 21 ) 


where the R and r are two principle radii of the curved 
droplet interface measured on the side of the non- 
spherical droplet, as shown Figure 13 Since the base 


of the droplet is almost flat and does not deform rapidly 
during the later stage of the drainage of the gas layer, 
the pressure in the droplet and in the gas layer can be 
regarded to be in equilibrium, and the pressure in the 
gas layer, can be approximated by the pressure in the 
droplet, pd, 


^'^deformation — 


p {vo COS 0) h 

P 


(25) 


The characteristic velocity is the normal velocity of 
droplet towards the plate, vq cos 6 >, while the charac¬ 
teristic length is the film thickness h as calculated 
from the Nusselt solution [43]. At an impact speed 
of 0.30 m/s, for instance, Wedeformation = 0.18, and 
Redeformation = 62, which indicate that the surface 
tension is very important in resisting the deformation of 
the film interface. Equations (24-25) also indicate that 


thicker films which correspond to high film Re are easy 
to deform. In addition, the relative importance between 
the viscous and capillary forces is quantified by 


Pg^l 




( 22 ) 


Cadeformation — 


Wedeformation COS 0 

R'Cdeformation 7 


(26) 


The lubrication force applied on the droplet by the liquid 
film via the gas layer can be estimated by multiplying the 
pressure in the gas layer and the base area of the droplet, 
Sb ^ 7r{R — r)‘^. Therefore, the lubrication force is 


- lubrication 


+ - ) 7r{R - rf 


(23) 


where the effect of the ambient pressure Pa in Eq. (22) 


is cancelled by the ambient pressure on the upper side 
of the droplet. For the bouncing event shown in Figure 
image analysis shows that the ratio of the lubrication 
force to the gravitational force ^lubrication/^gravity ^ 10- 
Therefore, the lubrication force by the gas layer on the 
droplet can be much larger than the gravitational force 
of the droplet, and can result in bouncing. 

At impact, both the gas-droplet and gas-film interfaces 
undergo deformation m whose severity depends on 
the system parameters. An estimate of the liquid film 
deformation is obtained by comparing inertial, capillary. 


which is in the range of 0.003-0.03 for different impact 
speeds (0.30-2.99 m/s), indicating that the capillary 
force is much stronger than the viscous force. 

During impact, at larger Re values of the liquid films, 
the average liquid films are thicker and more deformable 
since the restoring capillary force becomes progressively 
weaker. The deformation of the gas/liquid-film interface 
delays the drainage of the gas layer, preventing its 
rupture. In addition, the force exerted by the film on 
the droplet at high film Re, which is transmitted through 
the gas layer, pushes the droplet away from the impact 
region. These factors also could increase the probability 
of droplet bouncing. 

To further understand the collapse of the gas layer, the 
lifetime of the gas layer was studied, and a histogram 
of this lifetime is plotted in Figure The starting 
time was determined from the high-speed images, which 
show noticeable capillary waves on the surfaces of the 
liquid film upon apparent impact. The ending time was 
determined from the images which show the rupture of 
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FIG. 14. Histogram of the lifetime of the gas layer for three 
impact phenomena: bouncing, partial coalescence and total 
coalescence. The dimensionless parameters are Oh = 0.0021, 
We = 4, and Re in the range of 54-636. The droplet speed 
is 0.30 m/s, the flow rate of the falling him is in the range of 
0.9-10.7 1/min, and the droplet diameter is 3.3 mm. 


the gas layer indicated by the capillary waves on the 
droplet and on the liquid him, as explained in Figure 
For the bouncing phenomenon, the lift-off of the droplet 
away from the him was considered as the ending time. 
Due to this distinctive feature and the fast propagation 
of the capillary waves in the high-speed images, the life 
time of the gas layer can be determined at an accuracy 
of about 0.2-0.4 ms at the hlming frame rate 5000 
fps. As shown in Figure the bouncing and total 
coalescence phenomena correspond to long and short gas 
layer lifetimes, respectively. This indicates that delaying 
(accelerating) the rupture of the gas layer results in 
bouncing (total coalescence). This is because the speed of 
the liquid at the top region of the droplet can be retarded 
via droplet deformation during the thinning of the gas 
layer. Following rapid rupture of the gas layer, the fluid 
near the droplet apex has a higher chance of catching up 
with the fluid near the bottom of the droplet, preventing 
the formation of secondary droplets, and leading to total 
coalescence. 


III.8. Dimple formation and droplet ejection for 
large droplets 

The Ohnesorge number can be varied by changing 
the droplet size. We consider a case shown in Figure 
\TE\ characterised by Oh = 0.0016, We = 6.5, and 
Re = 54. These parameters are similar to those 
that were used to generate Figure which depicted 
a case of partial coalescence; the decrease in Oh and 
the slight increase in We were brought about by the 
increase in the droplet diameter from 3.3 mm to 5.2 mm. 
The case presented in Figure illustrates an impact 
accompanied by the droplet surface switching from a 
convex to a concave shape, followed by the ejection of 


a small droplet. Inspection of Figure [14^ shows that the 
droplet shape prior to impact is far from spherical, as 
expected for large droplet size, i.e., large Oh numbers. 
Upon apparent contact with the liquid film (see Figure 
[T^), capillary waves are formed and propagate along 
the surface of the droplet towards the top of the droplet 
and also on the surface of the liquid film (see Figure 
[T^-d). The liquid in the top region of the droplet is 
less affected by the impact than that near the bottom 
of the droplet and continues to move towards the film 
by inertia (see Figure [T^-f). Following rupture of the 
gas layer, capillary waves form and propagate on the 
surface of the droplet towards its apex. The interaction 
of these waves with the downwards-moving fluid from 
the apex leads to the formation of a hat-shaped droplet 
(see Figure p!^-h) followed by a dimple at the droplet 
apex (see Figure [T^-k). The formation of the dimple 
is facilitated by the relative weakness of capillary forces 
for such large droplets and their consequent inability 
to act as effective restoring forces to counterbalance 
severe interfacial deformation. For the concave interface 
associated with the dimple, capillarity pulls the liquid at 
the rim of the dimple causing it to shrink (see Figure 
p!^-n). The collision of the interface at the rim of the 
dimple results in the ejection of a small droplet (see 
Figure [T^-p). Due to the flow in the liquid film, the 
droplet is not symmetric, and the small droplet is ejected 
in a direction which is not perpendicular to the substrate 
but deviated in the streamwise direction. The liquid 
eventually settles down as the energy from the impact 
is dissipated gradually (see Figure [T^-r). 

III.9. The asymmetric splashing process 

We devote the concluding part of this section to quan¬ 
titatively examining the details of asymmetric splashing 
at different film flow rates, droplet speeds, and droplet 
sizes. The asymmetric splashing process was studied to 
consider the effects of relevant parameters and quantified 
by the propagation of the baseline of the droplet, as 
shown in Figure [T^. The streamwise direction of 
the falling film is defined as the positive direction of 
the movement of the baseline, and the initial contact 
point of the droplet on the falling liquid film is defined 
as the origin. The initial contact time is determined 
by interpolating the interface position from the image 
sequence, which provides a higher precision of the initial 
contact time than the time delay between two consecutive 
frames. The length of the baseline of the droplet, and the 
position of the front and rear point is normalized by the 
diameter of the droplet, and the time is normalized by 
the speed and the diameter of the droplet: 

X = x/d (27) 

t = tvo/d (28) 

Since the droplet impact in the present study is different 
from the perpendicular impact of droplets on quiescent 
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FIG. 15. Dimple formation and droplet ejection during the 
impact of a large droplet on a falling film. The dimensionless 
parameters are Oh = 0.0016, We = 6.5, and Re = 54. The 
droplet speed is 0.3 m/s, the flow rate of the falling film is 
0.9 1/min, and the droplet diameter is 5.2 mm. See Ref. [6^ . 
Movie 5. 


films, the propagation of the front and rear of the 
baselines are asymmetric, as shown in Figure The 
degree of asymmetry is influenced by the oblique impact 
of the droplet and the flow rate of the liquid film. Upon 
impact, the front point of the baseline moves along the 
flow direction, and the rear point of the baseline moves 
backwards. In addition, the front point moves at a higher 
speed than the rear point. 

The effect of the droplet size on the asymmetric 
splashing process is shown in Figure \TEjp . The curves 
of the baselines for different droplet sizes collapse, even 
though the dimensional data show clear difference at 
different droplet sizes (see the inset plot of Figure [T^). 
For smaller droplets, the baselines propagate much more 
slowly at both the front and the rear of the droplets, and 
the splashing process is much shorter than that for large 
droplets. This is because the mass and the momentum 
introduced into the film by small droplets are much less 
than those by large droplets. 

The propagation of the baseline for droplet impact 
at different droplet speeds is shown in Figure [T^. 
The propagation of the baseline appears to be weakly 
dependent on the droplet speed. The dimensional plot 
(see the inset plot of Figure [T^) shows that with 
increasing the droplet speed, both the front and the 
rear points move faster in the downstream direction. 
This is mainly due to the fact that droplets with higher 
speeds (larger momentum) can overcome the resistance of 
impact and move faster downstream than droplets with 



^ 0 1 2 3 4 5 '^0 2 4 6 8 10 12 


tVo/d tVo/d 

FIG. 16. Effects of the droplet size, the droplet speed, and 
the film flow rate on the propagation of the baseline during 
droplet splashing, (a) Schematic diagram of the definition of 
the coordinate system, (b) Effect of the droplet size. The 
droplet speed is 1.87 m/s, the film flow rate is 0.9 1/min, 
and correspondingly. Re = 54 and We in the range of 156- 
254. (c) Effect of the droplet speed. The film flow rate is 4.1 
1/min, the droplet diameter is 3.3 mm, and correspondingly. 
Re = 245 and Oh = 0.0021. (d) Effect of the film flow rate. 
The droplet speed is 2.67 m/s, and the droplet radius is 3.3 
mm, and correspondingly. Oh = 0.0021 and We = 316. The 
inset shows two images of crowns at different film flow rates. 


lower speeds. 

The effect of the film flow rate on the propagation of 
the baseline is shown in Figure [T^. At a small flow 
rate, the liquid film is thin, and the liquid sheet of the 
crown is ejected close to the substrate. In contrast, 
at a high flow rate, the liquid film is thicker, and the 
liquid sheet of the crown is ejected close to the direction 
perpendicular to the substrate, as shown in the inset 
images of Figure p!6| i. Similar phenomena also exist for 
droplet impact on quiescent liquid films [9]. Therefore, 
at a low film flow rate, the movement of the liquid sheet 
of the crown can promote the movement of the droplet 
baseline. With increasing the film flow rate, the film 
thickness increases, and the motion enhancing effect of 
the crown decrease. Consequently, the speeds of the front 
and rear points reduce. If the film flow rate increases 
further, the flow speed in the falling film also increases 
to the point where the film simply transports the droplet 
downstream so that both the front and rear points move 
in the streamwise direction. 


IV. CONCLUSION 

In this paper, the oblique impact of droplets on 
inclined falling liquid films was studied using high-speed 
imaging. The results of droplet impact on flowing liquid 
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films showed unique features different from the impact 
process on quiescent liquids. Under different impact 
conditions, various impact phenomena were observed: 
droplet bouncing, partial coalescence, total coalescence, 
and splashing. An impact regime map was generated 
for the impact of droplets on falling liquid films. 
Splashing phenomena occur at large Weber numbers. 
We, while other phenomena manifest themselves at low 
We. Low Reynolds numbers of the liquid films. Re, are 
associated with partial coalescence, while the probability 
of bouncing and total coalescence increases with Re. 
Dimple formation and droplet ejection were observed 
for the impact of large droplets. Splashing processes 
were quantified through the propagation of the droplet 
baselines during impact, and the effects of the droplet 
size, the droplet speed, and the flow rate of the falling 
liquid film were analysed. The analysis of the lubrication 
force during the impact process shows that a higher 
flow rate in the liquid film produces a larger lubrication 
force, slows down the drainage process, and increases the 
probability of droplet bouncing. 

The study mainly focused on the fluid dynamics of the 
droplet impact on falling liquid films. There are many 


open questions related to droplet impact that are worthy 
of investigation, such as the direct measurement of the 
flow field during the impact process, the evolution of 
the gas layer, and the effects of droplet impact on the 
flow in the liquid films. In various relevant applications, 
such as annular flow, spray cooling, and spray painting, 
the impact processes are much more complex. For 
example, the impact of swarms of droplets on liquid films 
involves the interaction of numerous impact processes; 
the flow of the gas phase may induce turbulence in 
the liquid films and droplet entrainment on the film 
surface. This investigation will deepen our understanding 
of the physics involved in the phenomena of droplet 
impact, leading to a marked improvement in processes 
and applications that rely on these phenomena. 
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